Spin-torque driven ferromagnetic resonance in a nonlinear regime 
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Spin- valve based nano junctions incorporating Co|Ni multilayers with perpendicular anisotropy 
were used to study spin-torque driven ferromagnetic resonance (ST-FMR) in a nonlinear regime. 
Perpendicular field swept resonance lines were measured under a large amplitude microwave cur- 
rent excitation, which produces a large angle precession of the Co|Ni layer magnetization. With 
increasing rf power the resonance lines broaden and become asymmetric, with their peak shifting 
to lower applied field. A nonhysteretic step jump in ST-FMR voltage signal was also observed at 
high powers. The results are analyzed in in terms of the foldover effect of a forced nonlinear oscil- 
lator and compared to macrospin simulations. The ST-FMR nonhysteretic step response may have 
applications in frequency and amplitude tunable nanoscale field sensors. 



Spin-torque driven ferromagnetic resonance (ST- 
FMR) has proven to be a powerful technique to in- 
vestigate magnetic anisotropy, damping and spin-torque 
interactions in nanometer scale metallic ferromagnets 
P, [2, H, 0]. The method involves injecting an rf cur- 
rent into a laterally confined giant magnetic resistance 
or magnetic tunnel junction device that drives one of the 
magnetic layers into resonance. Homodyne mixing of the 
rf current and device resistance oscillations results in a 
dc voltage across the device. 

In our earlier studies, ST-FMR was conducted on spin 
valve junctions that incorporate Co|Ni multilayers that 
have perpendicular magnetic anisotropy and a thick Co 
reference magnetic layer (HQ. The Co|Ni layer magnetic 
anisotropy was measured with rf currents in the sub- 
mA range, well below the devices' spin-transfer switching 
threshold of ^ 8 mA. In these studies the rf current was 
modulated on and off at a low frequency (~ kHz), and a 
lock-in amplifier was used to measure the mixing voltage, 
V. Measurements were conducted in a linear response 
regime in which the homodyne signal was proportional 
to the rf power, i.e. where the field swept resonance line, 
V(H)/Itf was independent of 7 r f. Here J r f is the ampli- 
tude of the rf current. The resonance lines were sym- 
metric and magnetization precession angles were small 
(typically around 2° and no more than 6°). 

In this Letter we present ST-FMR studies in a non- 
linear regime, where the amplitude of the rf current is 
much larger. The resonance field and the line width of 
the Co|Ni layer have been determined with increasing rf 
current amplitude and the rf current is not modulated, 
enabling an investigation of the effect of field sweep direc- 
tion on the resonance lineshape. We analyze these results 
with a foldover theory of nonlinear oscillators that pre- 
dicts a shift in resonance field and a hysteretic response, 
with a lineshape that depends on the field sweep direction 
(Fig. [2(b)). We also compare our results to macrospin 
simulations. 

Pillar junctions with a lateral dimension of 50x150 
nm 2 (Fig. [2(a)) were patterned on a silicon wafer us- 
ing a nanostencil process [7|. Junctions were deposited 
using metal evaporation with the layer structure || 1.5 
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FIG. 1: (a) Sample layer structure and ST-FMR circuit, (b) 
The foldover effect: precession amplitude vs field at fixed 
rf frequency for three different normalized rf power levels, 
(3 = 0.2pc, p c and 3f3 c (f3 c = 8/(3^3)). In the last case the 
curve is multivalued, to the left of the resonance there is a 
low and high amplitude branch, leading to field sweep direc- 
tion dependent hysteresis (red arrows). The dashed part of 
the curve is an unstable oscillation. Inset shows the magne- 
tization configuration. The applied field is in the z direction 
and the magnetization of the Ni|Co layer precesses about this 
direction at an angular frequency oo and with a transverse 
magnetization Mt. 



nm Cr| 100 nm Cu| 20 nm Pt| 10 nm Cu| [0.4 nm Co| 0.8 
nm Ni]x 3 | 10 nm Cu| 12 nm Co| 200 nm Cu||, where 
the thin (free) layer is composed of a Co|Ni multilayer 
and the thick (fixed) layer is pure Co. 

ST-FMR measurements were conducted with a mag- 
netic field applied nearly perpendicular (within ~ 2°) 
to the sample surface at room temperature. An rf cur- 
rent generated by a high frequency source is amplified by 
20 dB and is added to a dc current using a bias-T (Fig. 
[2(a)). Positive dc currents are defined such that electrons 
flow from the free layer to the fixed layer. The magne- 
toresistance is 2.1% and the pillar resistance is 1.53 ft in 
the parallel state (Fig. 1(b) of Ref. [5]). The inset of 
Fig. [2(b) shows the magnetization configuration. With 
an applied field greater than 0.3 Tesla, the Co|Ni multi- 
layer is magnetized along z, perpendicular to the plane of 
the layers. However, since the Co layer has a large mag- 
netization, hqM s ~ 1.5 T, it remains largely magnetized 
in plane, as represented by the vector rhp. rf currents 
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FIG. 2: (a) ST-FMR signal as a function of the applied field at 
different rf currents levels for / = 16 GHz. J r f was 1.4, 2.1, 3.3, 
4.6, 6.2, 7.8, and 9.0 mA respectively. Each adjacent curve 
is offset by 10 jjlV. The resonance at 7 r f=1.4 mA is amplified 
by a factor of 2 above the background of 37 /xV. Inset: dV/dl 
vs I of at /io#=0.675 T, which is the resonance field at /=16 
GHz in the linear regime, (b) Magnetic field where the step 
jump occurs vs 7 r f. (c) Vpeak/^rf vs 7 r f. The peak voltage is 
measured with respect to the baseline background voltage. 



cause the magnetization of the Co|Ni multilayer, M, to 
precess with a transverse amplitude M t at an angular 
frequency uo. 

ST-FMR was conducted with a constant rf power, and 
a dc voltmeter was used to measure the mixing voltage. 
The applied field was swept up and down at a rate 0.12 
T/min. Fig. [21(a) shows an example of such measure- 
ments. In this case, the rf frequency / is 16 GHz and 
/ r f ranges between 1.4 and 9.0 mA. Each curve is offset 
by 10 /iV for clarity. At the lowest rf amplitude (J r f=1.4 
mA), there is a 37 /iV non-resonant background, which 
we associate with nonlinearities in our measurement cir- 
cuit. This offset varies with the rf amplitude, which we 
attribute to other parasitic rectification processes in our 
measurement circuits, such as slightly non-ohmic con- 
tacts and junction heating effects. At 7 r f =1.4 mA we are 
in the linear response regime, the lineshape is symmetric, 
and the resonance field is 0.675 T. The critical current 
for spin-torque- induced dynamics at this field is ~8.5 mA 
(inset of Fig. [2(a)). The lineshape begins to become 
asymmetric at J r f =2.1 mA. Starting from J r f=3.3 mA a 
step jump of the mixing voltage occurs on the lower-field 
side of the resonance. The maximum ST-FMR voltage 
increases as 7 r f increases, and the step-jump field contin- 
ues to decrease. 

The shift of the resonance peak with increasing power 
can be understood in terms of forced oscillator in which 
the frequency of oscillator depends on the oscillation am- 
plitude, i.e. the system is a driven anharmonic oscillator. 
The micromagnetic energy of the free layer can be writ- 
ten as: 
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where m z = M z /M and M eff = M - 2K 1 /(fjL M ) is 
the effective easy-plane anisotropy. The first term is 
the Zeeman energy and the second term is written in 
terms of M e ff, which contains contributions due to shape 
anisotropy (Mo) and the layer perpendicular magnetic 
anisotropy (K\). The resonance condition in this ge- 
ometry (applied field along z) is given by [8[, 2irf = 
7/io(^res — M e fim z ), where 7 is the gyromagnetic ratio. 
For a fixed frequency the resonance field is given by: 
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As the precession amplitude m t = M t /M G increases, 
the longitudinal component of magnetization, m z ~ 
1 — ra 2 /2, decreases resulting in a shift of the resonance 
to lower fields (/ioMsff ~ 0.2 T for our Co|Ni multilayer), 
as observed in our experiments. 

For small rf current amplitudes the lineshape is 
Lorentzian: 



(H - Ho) 2 + (AH)* 



(3) 



with a half width at half maximum of AH = 0^/(7/^0), 
where a is a constant that depends on the spin-torque per 
unit current and the magnetization direction of the fixed 
magnetic layer, rhp. a is the Gilbert damping constant. 
With increasing precession amplitude one finds 0, [TO, EH, 
[l2[ to first order in ra 2 : 
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(e + /3m t 2 ) 2 + l' 



(4) 



E = -fi Mom z H + -/i M M eff mj, 



(1) 



Here m t = Tn t /m t max is the transverse component of 
magnetization normalized to the maximum precession 
amplitude, m 2 max = aI 2 f /(AH) 2 ; e = {H - H )/AH. 
The resonance peak position is given by: e = — f3 — 
— aM e ffItf/[2(AH) 3 ]. So the resonance peak position is 
predicted to shift with the rf power, I 2 { . 

Fig. DIb) shows a plot of m versus field for three 
different rf amplitudes. There is a critical value of 
(3 C = 8/ (3 a/3) ~ 1.54, which corresponds to the foldover 
threshold for the excitation amplitude 7 r f, for which the 
left side of the resonance curve has a vertical tangent. For 
sufficiently large current, (3 > f3 c , the resonance curve is 
multivalued, with two stable/metastable branches differ- 
ing in amplitude that overlap on the low field side of the 
resonance in Fig. [TJb). The dashed part of the curve in 
Fig. DJb) is an unstable oscillation. As a result hysteresis 
is predicted in field swept measurements, as indicated by 
the red arrows in the figure. Hysteresis begins when the 
shift of the resonance peak is greater than approximately 
the resonance linewidth (H res — Hq=/3 c AH). 

The experimental data (Fig. [2Ja)) have features qual- 
itatively consistent with the foldover model described 
above: (1) The resonance lines become asymmetric with 
increasing power; (2) There is a step jump in ST-FMR 
voltage above a critical rf amplitude; and (3) The reso- 
nance peak position decreases with increasing rf power 
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FIG. 3: Macrospin micromagnetic simulations of ST-FMR. 
(a) Voltage versus applied field for various rf currents, showing 
the high amplitude branch of the curve, (b) Resonance peak 
position and (c) V^eak/^rf vs 7 r f. The black squares refer to 
the low amplitude branch and the red circles refer to the high 
amplitude branch of the resonance curve. 

(Fig. [2Kb)). However, there are discrepancies with the 
model, particularly the fact that little hysteresis is ob- 
served at large rf amplitudes, when the resonance shift 
far exceeds the resonance linewidth, ji^AH = 0.02 T. 
Also, the resonance field shifts nearly linearly with in- 
creasing rf current. 

To further understand the data we have conducted 
macrospin simulations of the magnetization dynamics 
with an applied rf current. We take the micromagnetic 
energy given by Eq. [T] and an angular dependence of 
the spin-torque and magnetoresistance determined using 
continuous random matrix theory (CRMT) for our device 
layer stack [l3[. This model has been found to give spin- 
torques in agreement with those found in our experiments 
conducted in the linear response regime [5]. The torque 
is parameterized by tst = Ik(0)(h/2e)m x (m x m p ) 
with k(9) = (a + b sin 0)/(c + d cos 0) (a = 1.98, b = 
7.11, c = 28.74, d = 25.95). The results of our simula- 
tions are shown in Fig. [3j The high amplitude branch, 
with extends to lower field, is shown in Fig. [3Ka). The 
resonance lines shapes are similar to those observed in 
experiment. Fig. [3Kb) shows the shift of the peak posi- 
tion for both of high and low amplitude branches as well 
as the peak voltage. For large 7 r f the resonance peak po- 



sition shifts linearly with J r f and V pea ^/I r f saturates at 
large I r f (Fig. [31(c)), as seen in experiment. The main dis- 
crepancy is that larger rf currents are needed to produce 
the resonance shifts observed. 

It also appears that the resonance at high rf currents 
follows the high amplitude branch independently of the 
field sweep direction. The evidence for this is the fol- 
lowing. First, little to no hysteresis is observed. Second, 
the resonance shift far exceeds the linewidth as well as 
the resonance shift found for the low amplitude branch in 
micromagnetic simulations (Fig. [3Kb)). We suspect that 
there are thermally driven transitions between these dy- 
namical states that excite the high amplitude branch and 
reduce the hysteresis. The rf current may also heat the 
device reducing the magnetization of the free layer, in- 
creasing thermal fluctuations and leading to a larger shift 
of the resonance peak position with rf current than seen 
in the macrospin model. We also note that our analy- 
sis and macrospin simulations assume that the damping 
does not depend on precession amplitude. However, it 
has been noted that the damping may be nonlinear and 
this may play a role in understanding the large amplitude 
ST-FMR driven magnetization dynamics [1J]- 

It is interesting to estimate the maximum precession 
angle, # m ax- This can be done using the voltage peak 
as well as maximum resonance shift. The peak mixing 
voltage is given by V=jARI r f sin 77 sin# max , where AR 
is the junction magnetoresistance and rj is the angle be- 
tween the layers with no rf current applied (77 ~ 70°, in 
these experiments). The precession angle can also be es- 
timated from the shift of the resonance field using Eq. 
[21 Both approaches g lVe C/max ^ 

65°at J rf =9.3 mA, the 
largest rf current applied in the junction. 

In sum, these results illustrate that large rf currents 
can drive nonlinear magnetization dynamics, character- 
istic of any driven anharmonic oscillator. The observed 
nonhysteretic step response may prove useful for rf fre- 
quency and amplitude tunable nanometer scale field sen- 
sors. It will also be interesting to compare these re- 
sults to full micromagnetic simulations as well as to fur- 
ther explore the role of thermal fluctuations on ST-FMR 
foldover phenomena. 
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